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Abstract. The linear stability of the homogeneous equilibrium of non-relati- 
vistic fluids with mass flux and special relativistic fluids with the absolute 
value of the energy vector as internal energy is investigated. It is proved that 
the equilibrium is asymptotically stable in both cases due to purely thermo- 
dynamic restrictions; the only requirements are the thermodynamic stability 
and the nonnegativity of the transport coefficients. 



1. Introduction 

Thermodynamics is a theory of stabiHty. The existence of the entropy func- 
tion, a concave thermodynamic potential that increases in insulated systems is a 
stability requirement. Similarly, the existence of thermodynamic potentials with 
particular convexity properties and their tendency to a definite extremum in case 
of different homogeneous boundary conditions are stability requirements. This is 
a rather old, but essential point of view regarding the nonequilibrium thermody- 
namics of macroscopic systems, both continuous and homogeneous. The similarity 
of the mathematical structure of the continuum thermodynamic theories and the 
conditions in Lyapunov's direct method (see e.g. p]), was pointed out and inves- 
tigated by several authors |ll[3l[H[5]. Moreover, a complete and rigorous dynamic 
reinterpretation of ordinary equilibrium thermodynamics, the thermodynamics of 
homogeneous bodies from the point of view of stability is given in [6] (see also the 
references therein) . According to this point of view the homogeneous equihbrium 
of macroscopic continua, the thermodynamic equilibrium in case of homogeneous 
boundary conditions is asymptotically stable. Pure thermodynamic requirements 
are the only conditions without any further ado. In the following we call this 
essential property as generic stability. 

Let us give some clarifying remarks to the statement above. First of all generic 
stability can be violated in everyday situations. The loss of thermodynamic stabil- 
ity, that is the violation of concavity, the existence of phase boundaries, indicate 
that changes of the external conditions can trigger structural changes in the ma- 
terial. As a second remark let us emphasize, that generic stability is related to 
specific equilibria and not to steady states. The stability or instability of steady 
states is not related directly to generic stability [3j. Finally we emphasize that the 
existence of a Lyapunov function ensures that generic stability is nonlinear, and 
the linear stability of the same equilibria is only a necessary condition. 

This stability point of view of thermodynamics gives a sound and simple inter- 
pretation of the Second Law and shows clearly why the Second Law is fundamental 
for any theory of physics. Without this basic stability property of materials one 
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could not measure anything, repeated experiments could give rather different re- 
sults, small changes of the external conditions could drive physical systems far from 
their initial state. 

The problem of generic stability is well known and crucial for dissipative rela- 
tivistic fluids. The simplest relativistic generalization of the Fourier-Navier-Stokes 
system of equations, the theory of Eckart [7] is proved to be unstable and acausal 
and therefore not viable pj. The same has been proved for the so-called first order 
theories in general, where entropy four vector depends on the classical basic state 
variables [9]. 

The other conceptual issue of relativistic fluids, the requirement of causality, 
initiated a whole new family of dissipative relativistic fluids, where the basic state 
space is extended by different ways. In several cases the extension involves the dis- 
sipative part of the flux of energy-momentum as new variables of the theory. When 
the extension results in symmetric hyperbolic evolution equations, then causality 
is ensured. The price of the extension is the increased complexity of the equations 
and the larger number of necessary material parameters. There are a number of 
such kind of theories in the literature [101 El IH O El [15]. Only a few of them 
was investigated from the point of view of symmetric hyperbolicity and generic 
stability. 

Regarding the connection of causality and hyperbolicity there are arguments 
that from a physical point of view parabolic theories can be causal as long as 
they are stable [El El El E], because in that case the propagation speed of the 
vaHdity range of the theory can be smaller than the speed of light. On the other 
hand, stability and symmetric hyperboHcity are not necessarily connected, a stable 
parabolic theory can be extended. 

The example of the Israel-Stewart theory, the best known extended theory in 
the literature, demonstrates the rather involved relations of causality, stability and 
hyperbolicity of relativistic fluids. The original non-perturbed equations in these 
theories are not known to be symmetric (let alone causal) for arbitrary fluid states 
|20j . However, conditions of the symmetric hyperbolicity of the perturbation equa- 
tions are equivalent to the conditions of linear stability of the homogeneous equilib- 
rium. What is really disturbing here are these stability and hyperbolicity conditions 
themselves [211 122] (See Appendix A). These conditions give restrictions on the 
equations of state of the fluid beyond normal thermodynamic requirements, con- 
ditions for the newly introduced additional material parameters (e.g. they cannot 
be constant) and also restrictions regarding their relations. Moreover, the physi- 
cal content of these difficult relations is completely obscure. As the Israel-Stewart 
theory is the most popular in recent appHcations to the ffuid mechanic description 
of heavy-ion collisions all these conceptual problems are can have an experimental 
feedback [23, 24, 25, 26, 27, 24, 28, 29, 30j. 

Recently we have suggested a new approach to the problem. Our suggestion is 
based on a novel concept of internal energy [19]. The form of the internal energy 
was derived from a detailed investigation of the compatibility of the basic balances 
and the entropy inequality by modern thermodynamic methods [31]. We have 
investigated also the generic stability of the ffuid equations in the Eckart frame, 
and proved that thermodynamic stability and the positivity of the usual transport 
coefficients (heat conduction, viscosity) will ensure the stability without any further 
ado [T^. Moreover, we have given a simple extension of our internal energy that 
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gives relaxation equations for the dissipative fluxes similarly to the Israel-Stewart 
theory [32] . 

In this paper we investigate necessary conditions of the asymptotic stability 
of homogeneous equilibrium of insulated simple one component fluids. First we 
will show the linear stability of a generalized form of one component nonrelativistic 
fluids with a diffusion flux. These fluids were recently promoted by Brenner [33'j '34] 
and their thermodynamic suitability was shown in the frame of GENERIC 
and Classical Irreversible Thermodynamics [35]. The existence of diffusion flux 
in nonrelativistic fluids is equivalent to a general flow frame in relativistic fluids, 
when the velocity fleld of the continuum is not necessarily flxed to a conserved 
quantity. Then we investigate the generic stability of our relativistic dissipative 
fluid theory. We show that the conditions of the linear stability in a general frame 
are the thermodynamic stability and the positivity of the transport coefficients (heat 
conduction, viscosities and diffusion) without any further ado. Then we summarize 
and discuss our results. 

2. NON-RELATIVISTIC FLUIDS 

Here we develop a generalized Fourier-Navier-Stokes theory, from the point of 
view of non-equilibrium thermodynamics. Because we want to give a parallel treat- 
ment with the special relativistic case therefore we give the basic balances with the 
help of densities, introduce a conductive flux term in the balance of mass-particle 
number and apply a notation with indices in both cases. In case of non-relativistic 
fluids the indices are i, j, k, I, ... and are running from 1, ..,3. The dot denotes the 
substantial time derivative a = dta + v^dia, where dt is the partial time derivative. 

2.1. Basic equations. The balance of particle number will be given instead of the 
balance of mass 

(1) n + ndy + d,f = 0. 

In our one component fluid the particle number density n is related to the density 
as p = mn, where m is the molar mass. The balance of total energy is 

(2) e + ed^v' + dj' = 0. 

Here e is the conserved total energy density, P is the energy flux. The balance 
of the momentum, the Cauchy equation, will be written as 

(3) + p'djv' + dj {P'^ + fmv') = 0\ 

Here denotes the momentum density, we have assumed a continuum without 
internal moment of momentum, therefore P*^ , the pressure tensor is symmetric. 
There is an extra contribution to the momentum flux due to the diffusion flux. We 
assume that the momentum is related to the velocity and the mass as — pv"^. The 
internal energy density is the difference of the total and the kinetic energy densities 
as e = e — therefore the balance of internal energy is calculated from H])-© 
as 

(4) e + edy + a. (^l' - P'\ - = -P^'djv,. 

Here we can recognize the flux of the internal energy, the heat flow as 

2 

777 7; 

(5) q^^p^p^^^^^f_. 
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We assume local equilibrium, therefore the entropy is the function of the internal 
energy and the density. The Gibbs relation identifies the derivatives of the entropy 
and the potential relation defines the static pressure in our density based treatment 
(see e.g. 0). 

(6) de — Tds + ^dn, and e ~ Ts + — p. 

In the following we exclude phase transitions and assume thermodynamic stabil- 
ity, a concave entropy function. That can be ensured by the following inequalities 

^' deTdnT dnTdeT deT dnT \dnT ) " ' 

Here we exploited the equality of the mixed partial derivatives of the entropy 
^ = —■§^j^- Let us assume that the entropy fiux J* has the classical form [36] 

(8) J' = |-/f 

Let us remark that this form of the entropy flux is the consequence of the local 
equilibrium and flrst order nonlocality and can be derived from the Second Law 

m- 

Therefore the entropy balance can be written as: 



(9) s{e,n) + sdy + dir 



The usual thermodynamic fluxes are the heat flux and the diffusion flux and 
the dissipative pressure P'^ —pS^^ . The corresponding thermodynamic forces are the 
gradients of the entropic intensives di^ and — and the velocity gradient diVj 
respectively. Usually we assume a linear relationship between the thermodynamic 
fluxes and forces. For isotropic materials the representation theorems of isotropic 
functions give scalar coefficients in the Fourier-law, in the Fick-law and in the 
Newtonian pressure as 

(10) = Aa,l~xa.|, 

(11) f = xd.^-^d,^, 

(12) P'^ - pS'^ n*^' = -Tj{d'v^ + d^v' - ^OkV^'S'') - 'nAv'^S''. 

Here denotes the dissipative part of the pressure. A = \T^^ is the Fourier 
heat conduction coefficient, ^ is the diffusion coefficient and x is the coupling coeffi- 
cient of the Soret-Dufour effect, the coupHng between heat conduction and diffusion. 
The cross coefficients are equal according to the reciprocity relations of Onsager. 
77 is the shear viscosity and rj^ is the bulk viscosity, respectively. The inequality of 
the Second Law Q , that is the nonnegativity of the entropy production restricts 
the conductivity coefficients, therefore 

(13) A > 0, ^ > 0, - > 0, Tj>0, r]v> 0. 
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In the following investigations we will neglect the Soret-Dufour effect and assume 
that X = 0. 

2.2. Equilibrium. In thermodynamic (generic) equilibrium the substantial time 
derivatives and the dissipative fluxes g', and P*-' — p6^\ are zero. Therefore 
from (fTO| and ifTTj) it follows, that T(e,n) — const, and ^J.{e,n) — const.. The 
inequalities of thermodynamic stability ^ ensures the local invertibility of these 
functions, therefore in generic equilibrium e = const, and n = const., too. From 
(fT2| we get, that the velocity field is divergence and rotation free, that is djv^ = 
and diV^ — djV^ — 0. In the following we investigate the linear stability of the 
homogeneous equilibrium, where the relevant fields are constant in time and ho- 
mogeneous. Moreover, as the velocity field in our equations is a relative velocity 
related to an inertial observer, we may assume without restricting the generality 
that the equilibrium velocity field is zero. Therefore the equilibrium densities and 
fiuxes are 

n[x\t) = n = const., e{x-' ,t) ~1 ~ const., v^{x-' ,t) — , 

(14) f {x' , t) = 0^' , [x' , t) ^ 0^ , n'^' {x*' , t) = P'^ {x*' , t) - p{l, n)S'^ = 0'^ . 

2.3. Linearization. Let us denote the perturbed fields by {5n, Se, 5v^,5f,5q^, 5W^). 
The linearization of the balances H])-® around homogeneous equilibrium results 
in 

(15) = 5n + ndjSv^ + dj5j^ = 0, 

(16) = 6€+{-e+p)dj5v^ +dj5q\ 

(17) 0* = mnSv^ + d^p + d,6W^ = mn + ( ^d'Se + ^d'Sn] + d,SII'^ , 

\oe on J 

(18) 0^ = 6q^-Xd^^^5q^-X^-l-^d^Se+-^^d^Sny 

(19) 0^ = 5f + ^d^^ = Sf+^(^-^^^d'5e + ^^d^6ny 

(20) = +7^vdkSv''5'^ +7^{d'Sv' +d^v' + ldk6v''5'^). 

o 

The overline denotes the equilibrium value of the corresponding functions, e.g. 
p = p{e,n). 

In order to identify possible instabilities we select out exponential plane-wave 
solutions of the perturbation equations: SQ = Qqc^*'^'^^^ , where Qq is constant. As 
our equilibrium background state is a fiuid at rest, the substantial time derivative 
is the partial time derivative d/dt = dt- 

With these assumptions the perturbation equations follow as 

= T5n + ikn5v' +ik5f, 

= V5€+ {'e+p)ik5v'' + ikSq' , 

= VmnSv'' + ik ( ^de + ^Sn] + ik6U'='' , 
\de on J 

= TmnSyy + ikSIl'^y , 

= TmnJi;^ + 
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(21) 



= Sq"^ - ikXl ——5e+ ——5n 
^ \deT dnT 

= 5qy = 5q\ 



5jy = sf, 

SIF^ + ik + rj,^ 

SW'y + ikTjSv^ , 
SU'"' + ik7jSv\ 
+ ikr]y6vy , 
5W + ikr^ydv^, 

5wy. 



We can put the equations above into the following matrix form 

(22) M\5Q^ = 0. 

Here SQ^ represents the list of 17 fields which describe the perturbation of the 
Fourier-Fick-Navier-Stokes fluid: 

SQ ^ ((5n,5e,(5w^,(5g^,m^^,,5j^; 

6uy',5qy,sq^sjy,6f). 

Then the 17x17 matrix M^^ can be written in the block diagonal form 



(23) 











o\ 
























I 








V 



where I is an 5x5 identity matrix and the submatrices R„r and N^^ are defined as 
follows 



(24) 



(25) 





/ mnT 


ik 0\ 










= 1 ikr] 


1 , 










\ikr]y 


ij 








r 





ikmn 








ik\ 





r 


ik{e + p) 


ik 












rp 





ik 







-ikX^j, 





1 





















1 








ikfi 





1 


0/ 



where f] = 8ri/3 + rjv and we have removed the overline above the equilibrium 
quantities. 
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Exponentially growing plane-wave solutions of ([55| emerge whenever F and k 
satisfy the dispersion relation 



(26) 



detMnr = (detN„r)(detR„r)^ = 



with a positive real F. The roots of this equation are the roots obtained by setting 
the determinants of either N„r or rt^-p to zero. 
The determinant of R„r gives the condition 



pT + rik' ^ 0, 



which results in a negative F. 

The determinant of N„r gives the following dispersion relation 



■q-Xp 



d 1 



deT 



■ A^mnA 
A:''Aer?A + A:*^(e+p) 



9 .91 

dp d p, dp d fi 
'dednT ~ dndeT 



d 

+ (e+J5)^ + 



dp 



k^mnX 



de 
dp d 1 
'dednT 



F+ 



dn 
dp d 1 
' dndeT 



= 0. 



where 



deTdnT deTdnT - ' 



due to the concavity of the entropy. 

According to the Routh-Hurwitz criteria [38], the real parts of the roots of a 
third order polynomial a^x^ + aix^ + 022; + 03 = are negative whenever 



(27) 



ao > 0, 
ai > 0, 
aia2 — aoas > 0. 



We can see, that the first two conditions of fGOl are fulfilled, because the posi- 
tivity of the density and the transport coefficients and the concavity of the entropy. 
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The third condition after some calculations gives the following expression 



aia2 - aoas 



X d 1' 
rj + Xmn— — 
oe 1 



77-^ — ^ - 2ri\mn— — - (mn)'^A^ — — A 

' ^ ' dnT deT ^ ' "^'^ 



dnT 



4 mn 



d n 



77— — + mnAA 



d_l_ 

deT 

d 1 



d fi 



1 2 



dnT dnT 



9 ly d 1 

mn\ — — 77— 

^ deT 'deT 



d 1 

+ (e+p)2mn—- 



irnnY 



d ^ 



d 


r 


deT 


)] 


) 







> 



One can see, that the first four lines of the above expression are positive, term 
by term due to the thermodynamic conditions ([7]) and l|13p . The positivity of the 
last two lines can be verified if we recognize that the expression is a second order 
polynomial of (e + p). The coefficient of the quadratic term is positive and the 
discriminant of the polynomial simplifies to the following form: 

Ar]{mnf (-ry + mnA^;^ ) A < 0, 



deT 



As the discriminant is negative, the polynomial is positive. 

Therefore we have proved the linear stability of the thermodynamic equilibrium 
of Fick-Fourier-Navier-Stokes fluids. The linearized equations around the homoge- 
neous equilibrium of these fluids are asymptotically stable. 



3. Special relativistic fluids 

3.1. Balances of particle number, energy and momentum. For the metric 
(Lorentz form) we use the 5"^ = diag{—l, 1, 1, 1) convention and we use a unit speed 
of light c = 1, therefore for a four- velocity m° we have Uau"" = — 1. A°^ = g°'j^ + u°'Ub 
denotes the u-orthogonal projection. We perform the stability investigations in a 
general frame independently of the conventions used by Eckart or Landau and 
Lishitz. 

The particle number flow can be expressed by the local rest frame quantities as 

(28) = nu" + f. 

Here n = —UaN°- is the particle density and j° — A^^N^ is the particle flux in a 
comoving frame. 

The particle number conservation is described by 

(29) daN'' = h + ndaU'' + daf = 0, 

where n—^~ u^dan denotes the derivative of n with respect to the proper time 

T. 
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The energy-momentum density tensor is given with the help of the rest-frame 
quantities as 

(30) T'^'' = eu^u'' + u^g*" + + P°^ 

where e = UaUbT"'^ is the density of the energy^ = — UqA^^T"''' is the energy flux 
or heat flux, — — Uf,A°^T'^'' is the momentum density and P""^ = A^A^T'^'' is the 
pressure (stress) tensor. The momentum density, the energy flux and the pressure 
are spacehke in the comoving frame, therefore Ua9° = and UaP"^^ = UbP^'' = O''. 
Let us emphasize that this form of the energy-momentum tensor is completely 
general, but expressed by the local rest frame quantities. The energy-momentum 
tensor is symmetric, because we assume that the internal spin of the material is 
zero. In this case the heat flux and the momentum density are equal. However, 
the difference in their physical meaning is a key element of our train of thoughts. 
Heat is related to dissipation of energy but momentum density is not, therefore this 
difference should appear in the corresponding thermodynamic framework. 
Now the conservation of energy-momentum dbT"''' = is expanded to 

(31) dbT"'' eit" + eu'^dbu'' + eu" + u'^dbq'' + q^dbu" + + q'^dbu'' + dbP''^ = 0. 
Its timehke part in the local rest frame gives the balance of the energy e 

(32) - UadbT-''' - e + edaU"- + daq" + q^Ua + P'^'dbUa = 0. 

The spacelike part in the local rest frame describes the balance of the momentum 

(33) A^.^bT"'' = cm" + q'^dbU^ + q^dbu"" + A^.g'^ + A^^^P''" = 0. 

4. Thermodynamics 

The entropy density and flux can also be combined into a four-vector, using local 
rest frame quantities 

(34) 5" = -I- r, 

where s = —UaS'^ is the entropy density and J"" — S"" — u°'s — A'^j^S^ is the entropy 
flux. The entropy flux is w-spacelike, therefore UaJ"" = 0. Now the Second Law of 
thermodynamics is translated to the following inequality 

(35) 9a^" = S + sdaU" + daJ" > 0. 

Relativistic thermodynamic theories assume that the entropy is a function of 
the local rest frame quantities, because the thermodynamic relations reflect general 
properties of local material interactions. The most important assumption is that 
the entropy is a function of the local rest frame energy density, the time-timehke 
component of the energy momentum tensor according to the velocity fleld of the 
material [7]. Deflnitely the thermodynamics cannot be related to an external ob- 
server, therefore the dependence on the relative kinetic energy is excluded. This 
interpretation of e in l(30l) is supported by the form of the energy balance (|32l) , where 
the last term is analogous to the corresponding internal energy source (dissipated 
power) of the nonrelativistic theories. 

In nonrelativistic fluids the internal energy is the difference of the conserved 
total energy and the kinetic energy of the material. However, also in nonrelativistic 
theories the constitutive relations must be objective in the sense that they cannot 
depend on an external observer, the thermodynamic framework should produce 
frame independent material equations. (This apparent contradiction of classical 
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physics is eliminated by different sophisticated methods and lead to such important 
concepts as the configurational forces or/and virtual power [39l [40l |4l]). However, 
without distinguishing the energy related to the flow of the material from the total 
energy one mixes the dissipative and nondissipative effects. The wrong separation 
leads to generic instabilities of the corresponding theory. 

Our candidate of the relativistic internal energy is related to the energy vector 
defined by E"- — —UbT"''' = eu°- + q°- . The energy vector embraces both the total 
rest frame energy density and the rest frame momentum. Therefore its absolute 
value e = ||£'|| = ^/—EaE°- = \/e^ — QaQ"' seems to be a reasonable choice of the 
scalar internal energy. Its series expansion, when the energy density is larger than 
the momentum density is analogous to the corresponding nonrelativistic definition 

(36) e=y/e^-qaq''~e-^ + ... 

2e 

Thermodynamic calculations based on the Liu procedure support this assump- 
tion [H]. Let us emphasize, that our candidate of internal energy is not related 
to any external reference frame, only to the velocity field of the material. In a 
Landau-Lifschitz frame the energy vector is timehke. 

Assuming that the entropy density is the function of the internal energy and the 
particle number density s{e,q°',n) — s(-\/ — ga?", n) leads to a modified form of 
the thermodynamic Gibbs relation and the potential relation for the densities as 
follows 

(37) de — —dqa = Tds + iidn, and e — ^ ^° — Ts ^ p + /in. 

e e 

Here T is the temperature, p is the pressure and /i is the chemical potential. 
Equivalently the Gibbs relation gives the derivatives of the entropy density as fol- 
lows 



ds 
de 



1 ds 
T' dq^ 



{e,n) 



qa_ ds_ 

eT ' dn 



T' 



For the entropy flux we assume the classical form 

(38) r = ^l-t^. 



Now we substitute the energy balance l(32|) and the particle number balance l(29|) 
into the entropy balance l(35|l and we arrive at the following entropy production 
formula: 

das'" = s{e,q\n) + sdaU'' +dar 

- ^e + —r + —n + sdu'^ + d^~d^^ 
" de^ dq-^ ^ dn^ "^ "T '^"T 

= -^{edaU'' + daq'' + q'^Ua + P''^ dbUa) - ^qa + sdaU'' 

1 1 e 

+ ^{ndaU'^ + daf) + q'^da^ + ^daq'' - fda^ - ^d^f 

(39) = -/U^-^(p"^-(p+^)a"^)9.u, 
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According to this quadratic expression and the potential relation in l(37l) the 
viscous pressure is given by 



jab ryab I ^ ^ \ \ab 



Therefore the (|40|) entropy production can be written as 

(40) ~ fda^ ^n'^'daU, + (^54 - I - > 

In isotropic continua the above entropy production results in the following con- 
stitutive functions assuming a linear relationship between thermodynamic fluxes 
and forces 

(41) g'^ = AA'^''(^a,i-i-^ 

(42) r = -eA^^s.^, 

(43) n"^ = -2ry(9^i^) - 7/„ah^^^A°^ 

where the bracket denotes the symmetric traceless part of the spacelike tensor 

Here JS]), l(42|) and l(43|) are the relativistic generalizations of the Fourier law 
of heat conduction, the Fick law of diffusion and the Newtonian viscous pressure 
function. The shear and bulk viscosity coefficients, rj and 77^, the heat conduction 
coefficient A, and the diffusion coefficient ^ are non negative, according to the 
inequality of the entropy production l|4ip . We may introduce a relaxation time 
T = A/e in l(4T|) . as usual in second order theories. 

The equations l(29|) , l(32|) and l(33|) are the evolution equations of a relativistic heat 
conducting ideal fluid, together with the constitutive functions l(42|) . l(43|) and the 
relaxation type equation l(4T|) . As special cases we can get the relativistic Navier- 
Stokes equation substituting l(43|) into (|33|l and assuming 5" = 0, or the equations 
of relativistic heat conduction solving together (|4T|) and l(32|) assuming that IT*'' = 
and = const.. 



5. Linear stability 

In this section we investigate the linear stability of the homogeneous equilibrium 
of the equations (|29l) . l(32|) and l(33|) together with the constitutive relations l(4T|) - 
(|43| . Similar calculations are given by Hiscock and Lindblom both for Eckart fluids 
[5] and Israel-Stewart fluids [42] . 

5.1. Equilibrium. The equilibrium of the above set of equations is deflned by 
vanishing proper time derivatives and by zero entropy production with vanishing 
thermodynamic fluxes 



(44) 



n^*" = 0, f = 0, and q" = 0. 
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Therefore according to the balances and the constitutive functions the equihb- 
rium of the fluid is determined by 



(45) n = const, e — const. ^ T = const., ^ — const., p — const. 

(46) Sali" = 0, daUb + dbUa = 0. 



In addition to the above conditions we require a homogeneous equihbrium ve- 
locity field 



(47) Ua = const. 



5.2. Linearization. We denote the equilibrium fields by overline and the per- 
turbed fields by (5 as Q = Qq + SQ. Here Q stands for n, e, q"-, j° and 11"''. The 
linearized equations ([291), ^3^, (01]), |[42]) and |[43]) around the equiHbrium 

given by (|44l) - l(47|) become 



(48) = Sn + ndaSu" + daSf, 

(49) = Se+ie+p)dJu'' + daSq'', 

(50) = (e+p)(5u'^-h A"''9fc(5p + (5g« + A°9h(5^'=^ 

(51) ^ ^,"-AA-(^a.4-|^i-§), 

(52) = Sr+^A'^%S^, 

(53) = SW'' + T^^daSu'' A"'' + T^A'"'A'"^ {d^Sud + ddSu^-ldeSu'^Acd)- 



The perturbation variables satisfy the following properties inherited from the Hn- 
earization of the original ones 



= U°'Sqa = U°-5Ua = u'^SIVab = 5IVab - ^liba 



For the stability investigations we introduce exponential plane-wave solutions 
of the perturbation equations: 5Q = Qoe^*"'"*'^'^, where Qo is constant and t and 
X are two orthogonal coordinates in Minkowski spacetime. As our equiHbrium 
background state is a fluid at rest we put w"9a — dt. 
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With these assumptions the set of perturbation equations follow as 



(54) 



T6n + ikndu^ + ikSj^, 
T6e + (e + p)ik5u^ + ik5q^ , 

r(e + p)5w' + ik{dep6e + dnpSn) + T^g^ + ifcOT^^, 
r(e + p)6uy + rSqy + ikSW^, 

r(e + + rjg^ + 



1 



= Sq" - ikX de-Se + d„-6n + -TSu'' + —TSq=' 



A, 



A 



Te 



SqV + ^TSuy + ^T5qy, 
I 1 e 

5q^ + ^TSu^ + ^T5q\ 

Sf = Sf, 

kfjSu'^, 
krjSuy, 
krjSu^ , 
kfj^du^ , 



5W + 

+ 



Here we have introduced shortened notations for f] = r]y + ^r] and fjv = r}v — |r? 
and for the partial derivatives of the thermodynamic quantities as 9e = ^ and 
dn = -§^- We can put the equations above into the following matrix form 



(55) 



M\5Q^ = 0. 



Here 6Q^ represents the list of fields which describe the perturbation of the fluid: 

5Q = {5n, 5e, Sv", Sq", m^^, 5f; 
, Sqy , (5n^^ , SUyy Su\5q'',SU''\5W; 

Then the 17x17 matrix M can be written in the block diagonal form 



(56) 



M 



/N 0\ 

R 

R 

\0 1/ 



where I is a 3x3 identity matrix and the submatrices R and N are defined as follows 
(57) R = 



/{e + p)r 


r 


ik 


0\ 




1 + r— 








ikr] 





1 





\ ikfjy 








V 
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(58) 



N = 



/ 


r 





ikn 















r 


ik{e + p) 


ik 










ikdnP 


ikdeP 


T{e+p) 


r 


ik 









ikXde^ 


rA 


1 + — r 
















ikfj 





1 





V 




ik^dg ^ 











1/ 



Exponentially growing plane-wave solutions of (|55l) emerge whenever F and fc 
satisfy the dispersion relation 

(59) det M = (det N) (det R) 2 = 

with a positive real F. The roots of this equation are the roots obtained by setting 
the determinants of either N or R to zero. 
The determinant of R gives the condition 

.2 



eT 



P 



Te 



T + ffk^ = 



The real parts of the roots of this polynomial are negative because the coefScients 
of both the linear and the quadratic term are positive. 

The determinant of N gives the following dispersion relation 



Te 



p + k^ 



eT 



F2fc2 i^f, + (e + p)^d, 
Ffc^ ( (e + p)deP + ndnP + 



k^n 



\_ 

eT 
1 



'PCdn-, 

(ndnP 



pdeP + fc^77^9„ 



T 



k'^(d.^pdn^-d,^dnP 



eT 



't' 



Xr 



dePdnj; 



dnPde — 



+ ae+p) {depdnj; 



dnPdej; 



k^fjXS, 



de^dn — 



0. 



According to the Routh-Hurwitz criteria 



fourth order polynomial aox 



[38], the real parts of the roots of a 
- a^x + 04 = are negative whenever 



Go 
ai 

a4al 



> 
> 
> 
> 



0, 
0, 
0, 
0. 



aia2 

(60) (aia2 - 0003)02 

We can see, that the first two conditions of (|60l) are fulfilled according to the 
Second Law, the nonnegativity of the entropy production and the thermodynamic 
stability, the concavity of the entropy function. 

Let us observe that one could define two set of intensives according to the deriva- 
tives of the entropy by the total and the internal energies [32]. However, in the 
stability investigations we need only the equihbrium values of the thermodynamic 
relations, when g° — 0, therefore the two set of intensives as well as their partial 
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derivatives coincide. Hence the thermodynamic stabiHty is required in the same 
form as it was in the non-relativistic case, but the state variables are the particle 
density n and the total energy density e, instead of the density and the internal 
energy. ^ 



(61) 
(62) 

(63) 



de^ < 0, 

A := de^dnj; - dn^dej; < 0. 



As in the nonrelativistic case the following identities are useful 

deP = -{e + p)Tde^ + nTde^, 

dnP = - (e + p)Tdn ^ + nTdn j; , 

f) 1. - -r> H 
On rp ~ Oe rp- 

We can see that the first two inequalities of ([60l) are satisfied and the third 
inequality can be written as 



0102 — aooa 



?7(e+p) + [Xn"^ + S,{e + pf)dn 



e2T 



T 

ei] + X[ 11? - 2npdnTp, - np^d, 



1 



^X{2r^e{e+p)dn^ + Xep'{dn- 



T 

fc6 



+ 



(eT)2 "T 



Ap ( pdn- - nda- 



T 



This expression is positive term by term. Only the underHned part requires 
separate investigation. There one can recognise a second order polynomial of n 
with the negative discriminant 

4n2A < 0. 

In Appendix B we show that the positivity of the fourth inequality of ((60l) does 
not need any more condition beyond the concavity of the entropy and the positivity 
of the transport coefficients. 

Therefore we conclude that the homogeneous equihbrium of the relativistic heat 
conducting, viscous relativistic fluids is stable. We did not need to exploit any 
special additional stability conditions beyond the well known thermodynamic in- 
equalities. This is in strong contrast to the Eckart theory, where one encounter 
generic instabilities and to the Miiller-Israel-Stewart theory, where one should as- 
sume additional compHcated conditions f42j and Appendix A. Our statement is 
valid in arbitrary flow-frame, we did not need flxed the flow neither to the Eckart 
nor by the Landau-Lifshitz conditions. 
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6. Conclusion 

The suggested relativistic form of the internal energy depends on the momentum 
density, therefore the entropy function is the function of the momentum density, too. 
Moreover, in the first approximation we have a regular second order theory with only 
one additional quadratic term in the entropy four vector according to (|36|) . However, 
in contrast to any other extended theories there was no need to introduce additional 
parameters, the coefficient of the quadratic term, and therefore the relaxation time 
in the generalized Fourier equation, is fixed. Therefore our theory of dissipative 
relativistic fiuids can be considered as the minimal stable extension of the theory 
of Eckart. 

We have investigated the analogous non-relativistic theory to emphasize the dif- 
ferences and similarities of the conceptual questions. However, it is well known, 
that the suggestion of Brenner is untenable as it is, because the conservation of the 
moment of momentum requires the equality of the mass current and the momen- 
tum density, similarly to the relativistic case, where the time-spacelike and space- 
timelike parts of the energy momentum are equal in any frame as a consequence 
of the symmetry of the energy-momentum density in absence of internal moment 
of momentum 03l|44]. However, as the relativistic case demonstrates clearly, the 
question is not clarified completely, one can figure out several ideas to improve and 
further generalize non-relativistic fiuid dynamics e.g. as the nonrelativistic limit of 
the recent suggestion. 

There are arguments that a suitable chosen fiow-frame would eliminate the in- 
stabilities of the Eckart theory pi J.5J. However, the simplest choice to fix the fiow 
is the one of Eckart (coming from the concept of the non-relativistic baricentric 
velocity) with a vanishing particle fiux. The other natural choice, to fix the fiow to 
the energy fiux according to Landau and Lifshitz [46, 47\, leads to a heat fiux that 
is strictly connected to diffusion. Here we have shown that from a stability point 
of view the fiow frames are not necessary ingredients of a dissipative relativistic 
theory of fiuids. 
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8. Appendix A: Stability conditions of the Israel-Stewart theory 

The Israel-Stewart theory is based on the following form of the entropy four 
vector 



Here 11 — and tt"^ — (Jl°'j) are the trace and the symmetric traceless part of 
the dissipative pressure and /3o, /32 7 ^0: «i are the scalar valued Israel-Stewart 
coefficients. Stability conditions were calculated by Hiscock and Lindblom in [22j. 
The independent variables of the thermodynamic functions are different that we 
have used in the previous calculations. They use specific entropy s = s/n and find 
useful to introduce = ^ as independent variable. They require the positivity 
of the following quantities as conditions of the stability (formulas (53)-(60), with 
A = 1): 



(64) 






1 


de 




(e 


+ p) dp 


s 


(65) 






1 


de 


dp 


(e 


+ p) ds 




(66) 


^3 


= (e 


+ P) 


1 - 


de 
9P , 


(67) 


Qi 




+ P) 


2/32 + /5 


= (e 




2/3i/3i 


(68) 












(69) 


ne 


= f3i 


2 




2al 
3/32 


(70) 




= f3i 






(71) 




= P2 









2ai 



1 dT 



nT^ ds 



It is worth to express the partial derivatives of the thermodynamic quantities in 
the above expressions by partial derivatives with the energy density and particle 
number density as variables - using the previous shorthand notation - used in the 
our recent investigations: 



(72) 
(73) 

(74) 



de 
ds 



de 
dp 

dp 
ds 



dT 
'ds 



ie+p)T 



{e+p)^deT-n^T^dn§' 

,3 



[e+py 



n-T^ f n{e + p)dnT + n^Tdn^ 

''^\{e+prd,T-n^T^dn^ 
fj, — nTdn ^ 
{e+p)dnT + Tti-nT^dn^' 
T . 



-deT. 



n 
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With the graphic representation method of inequalities of the partial derivative 
of Hiscock and Lindblom [21] one can understand the physical content of the con- 
ditions. However, the meaning, the origin why one should require such a difficult 
set of conditions in addition to the thermodynamic restrictions, remains obscure. 



9. Appendix B: Fourth Routh-Hurwitz condition for relativistic 

FLUIDS 

The fourth inequality in l|60p can be written in the following form 



(aia2 - 0003)03 - 0^04 



{e+p)T] in^dnj; - 2n{e + p)dri^ - n{e + p)^de^ 



^ (-e{e+p)de^^ + in'^dn!^ - 2n{e + p)dn^ - n{e + pfde^ 



^[v\e+p)d^!^+ 
^ ^n2(3e + 2p)2(^9„i^ - nA{e + pf dnl;d,,^ + {e + pf {ie + p)2 [ d,. 



V 



{e+pfpde^dn^ 



ndn— - (e + p)dn— n 9„— - 2np9„— ~ p de— +epde — 



T 



T 



T 



/i\2 A 



e Vie+PY + -pie+p)d„^ (e +p)a„- - nd^^ 



T 



T 
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^) ( -ende^ + x(pde^+ndnl] 1 + 



eT 



+ 



3e(e + p) d, 



n dn- I (e-p) ( d„- ) - pde-dn- 



1 /i 



'^npdni; {{e+p)(dni;] + (e + 2p)de^dn^ 



p{e + p)de^ y{-e+p) (dn^j +pde^dn^ 
pdnj. ^2n2(5„^)'+2n(e + 2p)5„^9„|-(2e+p)(e+p) [dn^ + 



p{e+p)de^dn^ 



p'^ ( ndn^ - {e+p)dn 



I rp \ ' r / 11' fj-i 

e^p{e+p) (a^) 



rp p 

3 



er]d„- + Xp \{e-p){dn 



Pde^dnj; 



> 0. 



Here the first two underbraced expressions are identical second order polynomials 
of e+p with a discriminant 

Di = 4n^A < 0. 

Therefore the first part of the expression, the multiplier of /c^ is positive. Then 
the underHned expression is a second order polynomial of e, with a discriminant 

D2 = 4^ (ndn ^+pde^^ (^-{ndn^-pdn^y + 2p^A^ <0. 
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The second underlined and the third underbraced expressions are second order 
polynomials of n with negative discriminants 

D3=4{e+pf[dn^y(^p{Se + 2p)A-e{5e + Ap) j < 0, 

and 

L>4 = 4/A < 0. 

Therefore the coefficient of is positive, because the other terms and coefficients 
in that part of the expression are nonncgative. 

Similarly, in the coefficient of the three underlined expressions are second 
order polynomials of n with the following discriminants 



Da = 



4 (^dn^y (^2p{e + 2p)A - e(5e + 4p) (^5„^^ j < 0, 
4p (pA + e {dn^y^ (pA + ede^dn^^ < 0, 

Dr = 4 (dnff (^2p{e+p)A - e(3e + 4p) (^a„i^ j < 0, 

respectively. As every other term and coefficient is nonnegative in the expression, 
including the whole coefficient of A;^", we conclude that the whole expression is 
nonnegative, as we have already indicated at the end of the formula. 



[1 

p: 
[3: 
[4; 
[s: 
[6: 
[7: 

[9 

[lo; 

[11 

[12 

[13: 

[14 



References 

N. Rouche, P. Habets, and M. Laloy. Stability Theory by Liapunov's Direct Method. Springer 
Verlag, New York Inc., 1977. 

B. D. Coleman and V. J. Mizel. Existence of entropy as a consequence of asymptotic stability. 
Archive for Rational Mechanics and Analysis. 25:243-270, 1967. 

P. Glansdorff and I. Prigogine. Thermodynamic Theory of Structure, Stability and Fluctua- 
tions. Wiley-Interscience, London-etc, 1971. 

M. E. Gurtin. Thermodynamics and Stability. Archive for Rational Mechanics and Analysis, 
59:63-96, 1975. 

C. M. Dafermos. The Second Law of thermodynamics and stability. Archive for Rational 
Mechanics and Analysis, 70:167-179, 1979. 

T. Matolcsi. Ordinary thermodynamics. Akademiai Kiado (Publishing House of the Hungarian 
Academy of Sciences), Budapest, 2005. 

Carl Eckart. The thermodynamics of irreversible processes, HI. Relativistic theory of the 
simple nuid. Physical Review, 58:919-924, 1940. 

W. A. Iliscock and L. Lindblom. Generic instabilities in first-order dissipative relativistic fluid 
theories. Physical Review D, 31(4):725-733, 1985. 

W. A. Hiscock. Generic instabilities in first-order relativistic fluid theories. II. Havas-Swenson- 
type theories. Physical Review D, 86(6):1527-1532, 1986. 

I. Miiller. Toward relativistic thermodynamics. Archive for Rational Mechanics and Analysis, 
34(4):259-282, 1969. 

W. Israel. Nonstationary irreversible thermodynamics - causal relativistic theory. Annals of 
Physics, 100(1-2):310-331, 1980. 

W. Israel and J. M. Stewart. On transient relativistic thermodynamics and kinetic theory. II. 
Proceedings of Royal Society London A, 365:43-52, 1979. 

I-S. Liu, I. Miiller, and T. Ruggeri. Relativistic thermodynamics of gases. Annals of Physics, 
169:191-219, 1986. 

D. Pavon, D. Jou, and J. Casas- Vazquez. On a covariant formulation of dissipative phenom- 
ena. Annales de I'Institute Henry Poincare, Section A, 36(l):79-88, 1982. 



GENERIC STABILITY OF DISSIPATIVE NON-RELATIVISTIC AND RELATIVISTIC FLUIDai 



[15] H. C. Ottinger. Beyond equilibrium thermodynamics. Wiley-Interscience, 2005. 

[16] G. Fichera. Is the Fourier theorz of heat propagation paradoxical? Rediconti del Circolo 
Matematico di Palermo, XLI:5-28, 1992. 

]17] P. Kostdadt and M. Liu. On the causality and stability of the relativistic diffusion equation. 
Physical Reviews D, 62:023003, 2000. 

]18] V. A. Cimmelli. On the causality requirement for diffusive-hyperbolic systems in non- 
equilibrium thermodynamics. Journal of Non-Equilibrium Thermodynamics, 29(2):125-139, 
2004. 

Jl9] P. Van and T. S. Biro. Relativistic hydrodynamics - causality and stability. The European 
Physical Journal - Special Topics, 155:201-212, 2008. Zimanyi'75 Workshop Proceedings, 
larXiv:0704l0 39y2. 

]20] R. Geroch and L. Lindblom. Causal theories of dissipative relativistic fluids. Annals of 

Physics, 207:394-416, 1991. 
]21] W. A. Hiscock and L. Lindblom. Stability and causality in disssipative relativistic fluids. 

Annals of Physics, 151:466-496, 1983. 
[22] W. A. Hiscock and L. Lindblom. Stability in disssipative relativistic fluid theories. Contem- 
porary Mathematics, 71:181-220, 1988. 
]23] A. Muronga. Causal theories of dissipative relativistic fluid dynamics for nuclear collisions. 

Physical Review C, 69:0304903(16), 2004. 
J24] U. Heinz, Song H., and A. K. Chaudhuri. Dissipative hydrodinamics for viscous relativistic 

fluids. Physical Review C, 73:034904(11), 2006. 
J25] R. Baier and P. Romatschke. Causal viscous hydrodynamics for central heavy-ion collisions. 

European Physical Journal C, 51(3):677-687, 2007. [nucl-th/0610108' 
J26] G. S. Denicol, T. Kodama, T. Koide, and Ph. Mota. Shock propagation adn stability in causal 

dissipative hydrodynamics. Physical Review C, 78:034901, 2008. 
[27] G. S. Denicol, T. Kodama, T. Koide, and Ph. Mota. Stability and causality in relativistic 

dissipative hydrodynamics. Nuclear Physics G, 35(11):115102, 2008. 
[28] U. Heinz and H. C. Song. Causal relativistic hidrodynamics for viscous fluids. Journal of 

Phyics G, 35(10):104126, 2008. 
J29] Molnar D. and Huovinen P. Dissipative effects from transport and viscous hydrodynamics. 

Journal of Phyics G, 35(10):104125, 2008. 
J30] Molnar E. Comparing the first and second order theories of relativistic dissipative 

fluid dynamics using the 1+1 dimensional relativistic flux corrected transport algorithm. 

arXiv:0807.0544 

]31] P. Van. Internal energy in dissipative relativistic fluids. Journal of Mechanics of Materials 
and Structures, 3(6):1161-1169, 2008. Lecture held at TRECOP'07, arXiv:07121437 Jnucl-th]. 

]32] T. S. Biro, E. Molnar, and P. Van. A thermodynamic approach to the relaxation of viscosity 
and thermal conductivity. Physical Review C, 78:014909, 2008. arXiv:0805.1061 (nucl-th). 

]33] H. Brenner. Kinematics of volume transport. Physica A, 349:11-59, 2005. 

[34] H. Brenner. Navier-Stokes revisited. Physica A, 349:60-132, 2005. 

]35] D. Bedeaux, S. Kjelstrup, and H. C. Ottinger. On a possible difference between the barycnetric 
velocity and the velocity that gives translational momentum in fluids. Physica A, 371(2):177- 
187, 2006. 

J36] S. R. de Groot and P. Mazur. Non-equilibrium Thermodynamics. North-Holland Publishing 
Company, Amsterdam, 1962. 

J37] P. Van. Weakly nonlocal irreversible thermodynamics - the Ginzburg -Landau equation. Tech- 
nische Mechanik, 22(2):104-110, 2002. (cond-mat/0111307). 

J38] G. A. Korn and T. M. Korn. Mathematical Handbook for Scientists and Engineers: Def- 
initions, Theorems, and Formulas for Reference and Review. Dover, 2nd, revised edition, 
2000. 

]39] G. Maugin. The thermomechanics of nonlinear irreversible behaviors (An introduction). 

World Scientific, Singapure-New Jersey-London-Hong Kong, 1999. 
]40] M. E. Gurtin. Configurational forces as basic concepts of continuum physics. Springer, New 

York-etc, 2000. 

]41] M. Silhavy. The Mechanics and Thermodynamics of Continuous Media. Springer Verlag, 
Berlin-etc, 1997. 

]42] W. A. Hiscock and L. Lindblom. Linear plane waves in dissipative relativistic fluids. Physical 
Review D, 35(12):3723-3731, 1987. 



22 



PETER VAN 



[43] P. Kostadt and Mario Liu. Three ignored densities, frame/independent thermodynamics, and 
broken Galilean symmetry. Physical Review E, 58:5535, 1998. 

[44] Mario Liu. Comment on "weakly and strongly consistent formulations of irreversible pro- 
cesses". Physical Review Letters, 100(9):098901(1), 2008. 

[45] K. Tsumura and T. Kunihiro. Stable first-order particle-frame relativistic hydrodynamics for 
dissipative systems. 2008. arXiv:0709.3645 if2[nucl-th]. 

[46] L. D. Landau and E. M. Lifsic. Fluid mechanics. Pergamon Press, London, 1959. 

[47] Csernai L. Introduction to relativistic heavy ion physics. John Wiley and Sons, Chicester-etc, 
1994. 

Department of Theoretical Physics, KFKI, Research Institute of Particle and 
Nuclear Physics, Budapest, Hungary, and Department of Energy Engineering, Bu- 
dapest University of Technology and Economics 

E-mail address: vpethin8rmki.kfki.hu 



